Hadron Mass Spectrum from Lattice QCD 
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Finite temperature lattice simulations of quantum chromodynamics (QCD) are sensitive to the 
hadronic mass spectrum for temperatures below the "critical" temperature T c re 160 MeV. We show 
that a recent precision determination of the QCD trace anomaly shows evidence for the existence of 
a large number of hadron states beyond those known from experiment. The lattice results are well 
represented by an exponentially growing hadron mass spectrum up to a temperature T = 155 MeV. 
Using simple parametrizations we show how one may estimate the total spectral weight in these yet 
undermined states. 



Experimental data of multiparticle production in 
proton-proton collisions led Hagedorn to propose that 
the spectrum of hadronic states grows exponentially with 
mass [l|, Q. Such a spectrum arises naturally in the 
dual resonance model [3[ and, more generally, in mod- 
els of quark confinement, such as string models or bag 
models The spectrum of experimentally established 
hadronic states [5J is compatible with such an exponen- 
tial mass spectrum up to masses of approximately 1.7 
GeV Q . Most higher- mass hadron states are difficult to 
identify experimentally because of their increasingly large 
width and complicated decay properties. Nevertheless, 
significant efforts are being made to extend the baryon 
mass spectrum to higher masses Q , and searches for new 
meson states, including exotic states beyond those pre- 
dicted by the constituent quark model, are planned at the 
upgraded Jefferson Laboratory 12 GeV beam facility Q. 

A first-principles prediction of the hadron mass spec- 
trum from quantum chromodynamics by means of lattice 
QCD simulations would be highly desirable. Unfortu- 
nately, lattice simulations can only determine the masses 
of hadronic ground states and low excited states for given 
spin and parity Q . Here we point out that recent lattice 
simulations of QCD at finite temperature in the range 
of temperatures T < T c re 160 MeV [l(| are sensitive to 
the hadronic mass spectrum and permit to determine or, 
at least, constrain it beyond the experimentally estab- 
lished range. Earlier similar studies [11] which, however, 
were either based on lattice simulations with unphysically 
high quark masses and focused on the baryon sector of 
the hadron mass spectrum, or were performed for lat- 
tice actions that resulted in significantly higher values 
of T c [l2j . Other recent studies of high-mass resonance 
states were mainly aimed at their contributions to trans- 
port properties of QCD matter below the deconfinement 
temperature fl3l - [l5| . 

The so-called interaction measure (the QCD trace 
anomaly) , 



I(T) = (e-3P)/T 4 , 



(1) 



states. To see why this is so, we calculate I(T) for a ther- 
mal gas of non-interacting hadrons with mass spectrum 
p(m): 
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where we assumed that the hadrons are on mass shell: 
E 2 = p 2 +m 2 . Because the trace anomaly vanishes in the 
conformal limit, light hadrons contribute little to I(T). 
Heavy hadrons, on the other hand, contribute dispropor- 
tionally due to the presence of the factor m 3 in the inte- 
grand in the last line of eq. ([2]). In order to explore this 
quantitatively, we plot the integrand as a function of m 
for several fixed temperatures in the range T — 130 — 160 
MeV for an exponential mass spectrum of the form 



pirn) = cbe 



bin 



(3) 



derived from the trace of the stress-energy tensor = 
e — 3P, is especially sensitive to the presence of massive 



with b = (252 MeV)" 1 and c = 0.715. For T = 160 MeV 
(top curve) the integrand explores hadron masses much 
higher than for T = 130 MeV (bottom curve). Since the 
mass range of well established hadron states only reaches 
up to approximately 1.5 GeV for non-strange mesons 
and 2 GeV for baryons, one expects that the interaction 
measure I(T) is increasingly sensitive to experimentally 
unknown hadron states as the temperature exceeds 140 
MeV. 

This expectation is borne out by our results for the 
temperature dependence of the interaction measure. Fig- 
ure [H shows the results of a recent lattice calculation 
fioj as black dots (with error bars). The blue squares 
show, in comparison, the prediction of the hadron res- 
onance gas (HRG) model ([2} including the experimen- 
tally established states (23 non-strange mesons, 7 strange 
mesons, 19 non-strange baryons, and 17 baryons with 
open strangeness) , combining to a total of 795 states [E| . 
The red diamond symbols (with error bars) represent the 
difference between the lattice results and the HRG pre- 
diction. One notices that the HRG with the experimen- 
tally known mass spectrum provides an excellent descrip- 



2 




0.0 0.5 1.0 1.5 2.0 2.5 

m (GeV) 



3.0 



FIG. 1: (Color online) The integrand of eq. ([2]) in units 
of (1/GeV). The curves from bottom to top are for T = 
130, 140, 150, 160 MeV. Higher temperatures are seen to probe 
the hadron spectrum for increasingly larger hadron masses. 
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tion of the lattice results for temperatures up to T w 130 
MeV, but then falls below the lattice results. 

In order to understand the deviation we show, as the 
dashed (blue) curve, a prediction of the exponential mass 
spectrum with an upper cut-off at m c = 1.7 GeV. The 
curve is seen to be in excellent agreement with the pre- 
dictions of the HRG including only the experimentally 
known hadrons. The cut-off m c represents a reasonable 
compromise between the upper limits of the known me- 
son and baryon spectra. When we remove the cut-off and 
integrate without a limit placed on the allowed hadronic 
mass, we obtain the solid (black) curve, which follows the 
lattice results for the interaction measure up to T = 158 
MeV. The dotted (red) curve, which includes only masses 
above 1.7 GeV, provides a good representation for the dif- 
ference between the lattice results and the HRG curve. 
Again, this difference only seems to be describable by a 
nonintcracting hadron gas model up to a temperature of 
T = 158 MeV. 

The deviation of the HRG curve with an unlimited 
exponential mass spectrum from the lattice results for 
T > 155 MeV is not surprising, because this temperature 
approaches the quasi-critical temperature T c « 160 — 165 
MeV where deconfincmcnt and chiral symmetry restora- 
tion sets in [l6j . The HRG cannot provide a description 
of these phenomena. 

In order to explore the sensitivity of the interaction 
measure to the mass cut-off in the hadron mass spec- 
trum, we show in Fig. [3] the HRG predictions for the 
exponentially growing mass spectrum ^ with different 
cut-offs. The solid (black) curve shows the prediction of 
an unlimited spectrum. The dashed (blue) curve shows 
the prediction for an upper cut-off of m c = 1.7 GeV. The 
dotted curves show the results for upper mass cut-offs 
m c of 2.0 (bottom, blue), 2.5 (middle, green), and 3.0 
GeV (top, red), respectively. From the comparison of 



FIG. 2: (Color online) QCD interaction measure (e-3P)/T 4 
as a function of temperature T. The black dots (with error 
bars) show the results of a recent lattice calculation [IFJ . The 
blue squares show the prediction of the hadron resonance gas 
model including all experimentally established states. The red 
diamonds represent the difference between the lattice results 
and the hadron resonance gas prediction. The lines show 
the predictions from eq. (fj} for an exponential hadron mass 
spectrum. The solid (black) curve includes all masses; the 
dashed (blue) curve only masses up to 1.7 GeV; and the dotted 
(red) curve only masses above 1.7 GeV. The parameters of the 
mass spectrum are given in the text. 

these curves with the lattice results (black dots with er- 
ror bars), we see that at least a cut-off of m c = 2.5 GeV is 
needed to fit the lattice results up to T = 158 MeV. It is 
also clear that a significantly higher precision of the lat- 
tice data would permit to probe the cut-off dependence 
to even higher values of m c . However, Fig. [3] also clearly 
indicates that there is practically no sensitivity to masses 
above 3 GeV. 

To explore the sensitivity of our results to the analytic 
form of the resonance gas spectrum, we have repeated 
our calculation for a mass spectrum of the form originally 
proposed by Hagedorn [l| 

£ e rn/T H 

Ph( TO ) = T^~, 2T574- ( 4 ) 

The parameters were fit by comparing with the hadron 
resonance gas model prediction. An equally good fit as 
that in Fig. [2] is obtained for the parameters Th = 180 
MeV, A = 0.63 GeV 3/2 and m = 0.5 GeV. We will refer 
to this spectrum as the Hagedorn mass spectrum. The 
one difference is that in this case the cut-off mass must 
be chosen as m c = 1.9 GeV. Calculating the interaction 
measure with an unbounded distribution, i. e., setting 
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FIG. 3: (Color online) QCD interaction measure (e-3P)/T 4 
as function of temperature T. The black dots (with error bars) 
show the results of a recent lattice calculation [lj}. The lines 
show the predictions from eq. ([2} for an exponential hadron 
mass spectrum with different cut-off masses m c . The solid 
(black) curve includes all masses; the dashed (blue) curve only 
masses up to 1.7 GeV. The dotted curves show the results for 
mass cut-offs m c of 2.0 (bottom, blue), 2.5 (middle, green), 
and 3.0 GeV (top, red), respectively. 

m c = go, we obtain an identical prediction as in the 
case of Fig. [3]for the temperature (T « 155 MeV) where 
the Hagedorn mass spectrum begins to deviate from the 
lattice data. 

It may be surprising how two different parametriza- 
tions for the hadron mass spectrum can yield the same 
result for the temperature dependence of the interaction 
measure. This is due to two reasons: As pointed out 
above, in the relevant temperature range (T < 155 MeV) 
there seems to be very little sensitivity to resonances 
heavier than 3 GeV. Although the limiting temperature 
implied by the Hagedorn mass spectrum is considerably 
smaller (180 MeV compared with 252 MeV) than that 
associated with the exponential spectrum ([3]), both pre- 
dict nearly identical hadron mass distributions below 3 
GeV, as shown in Fig. 2) As a result, the predictions for 
the interaction measure from these two parametrizations 
differ only at very high temperatures, where they are no 
longer relevant to the lattice data. 

We conclude that recent lattice simulations of QCD 
with 2 + 1 dynamical quark flavors and physical quark 
masses already provide evidence for the existence of 
hadron states beyond those experimentally known. The 
interaction measure (JlJ, which is especially sensitive to 
high-mass hadron states, is in good agreement with an 
exponential mass spectrum of hadrons up to T = 158 
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FIG. 4: Number of hadron states below mass m shown on 
a logarithmic scale. The solid curve is the prediction of the 
exponential mass spectrum Q; the dashed line represents the 
prediction of the spectrum Q originally proposed by Hage- 
dorn. 

MeV, but a mass spectrum of the Hagedorn type can 
also provide a good fit. Increasingly precise future lattice 
simulations will enable a precision determination of the 
hadron mass spectrum beyond the range where individ- 
ual hadron states can be experimentally resolved. How- 
ever, our analysis demonstrates that the interaction mea- 
sure at temperatures below T c is not sensitive to hadron 
masses above 3 GeV. 

To obtain additional information about the hadronic 
mass spectrum, we propose the study of the variation 
of the interaction measure with baryon chemical poten- 
tial, strangeness chemical potential and/or quark flavor 
chemical potentials (l7j . A study in this vein was car- 
ried out in Ref. [HI, who used a hybrid model, i. e., a 
combination of known resonances below a mass cut-off 
and a Hagedorn spectrum above the cut-off, to explore 
the sensitivity of the energy density and the pressure to 
the parametrization chosen. However, to the best of our 
knowledge no study of the dependence of the interac- 
tion measure on quark chemical potentials has been been 
carried out. While we do not present any quantitative 
predictions in this Letter due to the scarcity of lattice 
data on this topic, the results of such a study may be 
easily discerned: The introduction of a baryon chemical 
potential will require one to decompose the hadron dis- 
tributions in Eq. J2]) into a mesonic and baryonic part 
with the Boltzmann distribution of the baryonic piece 
modified to include the baryon chemical potential, i. e., 
cxp(— E/T) — > cxp(— [E — fj]/T) with the implied oppo- 
site sign for anti-baryons. Derivatives of the interaction 
measure with respect to /i will allow for an estimation 
of the portion of these unknown states which lie in the 
baryon spectrum. Even more detailed analyses may be 
performed by the introduction of separate chemical po- 
tentials for each quark flavor. 
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